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Abstract. We present the latest generation of superconductor-insulator-fcrromagnct- 
superconductor Josephson tunnel junctions with a step-like thickness of the ferromag- 
netic (F) layer. The F-layer thicknesses d\ and d,2 in both halves were varied to obtain 
different combinations of positive and negative critical current densities j Ct i and j C) 2- 
The measured dependences of the critical current on applied magnetic field can be well 
described by a model which takes into account different critical current densities (ob- 
tained from reference junctions) and different net magnetization of the multidomain 
ferromagnetic layer in both halves. 



PACS numbers: 74.25.Fy 74.45. +c 74.50. +r, 74.70.cn 
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Superconducting spintronic elements, made up by superconducting (S) and 
ferromagnetic (F) layers, may improve future classical and quantum computing devices 
[1]. The underlying physics is based on the difference in spin-order (antiparallel 
alignment in S and parallel in F): the Cooper pair singlet injected into a ferromagnet 
gains a finite center of mass momentum, which, in turn, leads to an oscillating phase of 
the superconducting order-parameter [2]. In the ground state, the conventional, e.g. SIS 
(I: tunnel barrier), Josephson junctions (JJs) have the same sign of the order parameter 
in both superconducting electrodes, whereas in SFS or SIFS JJs the order parameters 
may have opposite signs [3, 4, 5, 6, 7], i.e., one deals with a u ir JJ", as the Josephson 
phase is <p = ir in the ground state due to a negative critical current J c < 0. The 
conventional JJs with I c > are called "0 JJ" in this context, as they have = 
in the ground state. Generally speaking I c is an oscillating and decaying function of 
the F-layer thickness dp, see Fig. 1(a). By using an F-layer with a step-like thickness, 
d\ in one half and d 2 in the other half, as shown in Fig. 1(b), one may fabricate a J J 
with different J Cj i = I c (d\) in one part of the JJ and = Ic(d 2 ) in the other part. 
If 7 Cj i and I c> 2 have different signs, one obtains a 0-7r J J [S], in which and 7r ground 
states compete with each other. At certain conditions [9] the ground state of the JJ is 
double degenerate, corresponding to a vortex of supercurrent circulating either clock- 
or counterclockwise, and thereby creating a spontaneous magnetic flux [10] |$| < $o/2, 
where $o is a magnetic flux quantum. 

Various 0-0 (/ c ,i,/ c , 2 > 0), vr-vr (I C)1 ,I Ct2 < 0) and 0-vr (7 Cjl > 0,7 Cj2 < 0) SIFS JJs 
based on diluted NiCu alloy were already fabricated and studied by us [8, 11, 12, 13]. 
In comparison with other types of 0-n JJs, such as SFS or d-wave/s-wave [14, 15], the 
intrinsic capacitance of SIFS JJs makes these JJs underdamped, so that one can study 
the dynamics, make spectroscopy [16] or use them in macroscopic quantum circuits [1]. 

In this paper we present a new generation of SIFS (Nb|A10a;|NiCu|Nb) JJs with a 
step in F-layer thickness. In comparison with the previous process the new one provides 
more superior JJs. First, additional Al interlayers within the bottom Nb electrode 
decrease the SI interface roughness. Second, computer control of the F-layer step etch 
enhances reproducibility. Third, a lower F-layer gradient gives a better sample yield 
[17]. 

The multilayer is deposited by DC magnetron sputtering on thermally oxidized 4-inch 
Si substrates. The 160 nm thick Nb bottom electrode, made up by four 40 nm Nb 
layers, each separated by 2.4 nm Al layers to reduce roughness, was covered by a 5 nm 
thick Al layer and thermally oxidized for 30 min at 10~ 2 mbar residual oxygen pressure 
and room temperature. To obtain many structures with different F-layer thicknesses 
in one fabrication run, we deposit a wedge-shaped F-layer (i.e. NiCu) alloy in order to 
minimize inevitable run-to-run variations. The F-layer gradient was reduced by a factor 
of ~ 2 compared to the previous process [11]. The multilayer stack was covered with a 
40 nm Nb cap layer. 

The stepped junctions were patterned using a four-level optical photolithographic mask 
procedure including SF 6 reactive etching and Ar ion-beam milling. The junctions were 
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partly protected with photoresist to define the step location in the F-layer, followed 
by i) selective Nb removal via SF 6 reactive etching, ii) Ar ion-etching of the NiCu and 
iii) in-situ deposition of Nb. The insulating layer between top and bottom electrode is 
self-aligned by ion-beam etching down to the AlO^ tunnel barrier and anodic oxidation 
of the bottom Nb electrode. Finally the top wiring is deposited. 

Various junctions were placed on the wafer within a narrow row perpendicular to the 
gradient in the F-layer thickness and were replicated along this gradient. One row 
contained a triplet of junctions including: 

• reference JJ with etched, uniform F-layer of thickness d\ 

• reference ir JJ with as deposited, uniform F-layer of thickness d 2 

• stepped 0-7T J J with step Adp in the F-layer thickness from d\ to d 2 

Our data are compared with a simple model, which takes into account different 
Ic,i 7^ Ic,2 and different net magnetizations Mi ^ M 2 in each half. Both the multidomain 
structure of the F-layer, yielding a stochastically distributed local magnetization, and 
the difference in magnetic thickness in both halves result in Mi 7^ M2. 

We consider a SIFS J J, shown in Fig. 1(b) with J C) i 7^ J C)2 . The fluxes $m,i 7^ $m,2, 
created by in-plane F-layer net magnetizations Mi and M2, are added to the flux $ 
provided by an external uniform magnetic field H . Note that Mi and M2 can be 
independently aligned either parallel or antiparallel to the applied field. Thus, their 
sign can be either positive (parallel alignment) or negative (antiparallel alignment). 
The restriction to ID magnetization is experimentally justified for weak magnets, the 
more general case of 2D in-plane magnetization in SIFS JJs is discussed elsewhere [18]. 
In the following we use an index i — 1, 2 to refer to two different parts of our J J. The 
magnetic flux density is fej = fi H2X L + $ Mii /Lj with the London penetration depth A^ 
and the length Li of each part. 

For a short JJ of length L = Li + L 2 and width w, both < 4Aj (Aj is the Josephson 
penetration length), the local phase is <pi(x) = <f)o + 27r/$o ■ hx with an arbitrary initial 
phase 0o and the maximum Josephson supercurrent for each H is given by 

I C {H) = max 

00 

where = j c ^wL. Further we assume that L\ = L 2 = L/2. 

The effect of the input parameters, such as I c ^ and 5>M,i, on I C (H) was discussed 
for our first 0-n junction (J c i > 0, I c2 < 0) in Ref. [13]. Now, using our latest 
generation of samples, we check the applicability of this model in the whole range of 

W/ C ,i = -1 . . . + 1. 

Fig. 1(a) shows the experimentally measured \I c (dF)\ dependence for reference JJs 
without an F-layer step. For JJs with an as-deposited (non etched) F-layer the I c (dp) 
vanishes at dp ~ 6.8 nm due to crossover from to it ground state, i.e. the sign change 
of 7 c ((ii?). The etched away F-layer thickness AdF ~ 1.7 nm was determined by the shift 
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of |/ c (djr)| dependence along dp axis after etching. Note that we have two reference JJs 
for each dp- the one was etched (dp = d\) and another left as deposited (dp = d 2 ). 

For J J with a step in the F-layer, Eq. 1 predicts the I C (H) patterns shown as color 
plot in Fig. 2 for I c ,2/I c ,i = —1 ... 1. The full range I C) 2/I c ,i = — 00...00 is included by an 
index change. The fluxes $m,i were set to zero. Experimentally the I C (H) dependences 
are measured at 4.2 K for three JJs with I c ,2/Ic,i < ( 0-tt JJs) and for four JJs with 
I c ^/I c ,i > (0-0 or Ti-71 JJs). They were taken on zero- field cooled samples for both 
field sweeping direction. The magnetic field was applied in-plane and parallel to the 
F-layer step, see Fig. 1 (b). The measured curves are shifted in Fig. 2 according to 
their I Cj 2/I c ,i ratio. At zero magnetic field the I c (0) changes from a global maxima at 
I c ,2/Ic,i = 1 to a global minimum at I c ,2/Ic,i — ~ 1 ; as the integrals in Eq. 1 gradually 
cancel out. For applied magnetic flux equal to the odd multiples of $0 the supercurrents 
cancel out completely only for the JJ with uniform I c (I Ct2 /I c ,i = 1)- For even multiples 
of $0 the i c = results from cancelation of supercurrents separately in each part of the 
JJ, irrespectively of the I c ,2/h,i ratio. Thus, Fig. 2 confirms the qualitative agreement 
of experimental data with the theoretical model taking into account l c \ 7^ l c2 . 

In Fig. 3 we show that a good quantitative agreement can be obtained by also 
allowing for independent magnetizations in both halves, i.e., $m,i 7^ $m,2- Although 
our analytical model can be applied to 0-0, tt-tt or 0-tt JJs, we focus on the 0-tt JJs 
as they are more sensitive to parameters and more interesting for applications. Three 
sets of JJs have the F-layer thicknesses d 2 = 7.3, 7.68 and 8.13 nm. Their lateral sizes 
(100 x 50 /iin 2 ) were comparable or smaller than Aj. The longest sample (0 JJ of set 
d 2 = 7.3 nm) has Aj ~ 75 fim. The measured I C (H) of the 0-tt and the and tt reference 
JJs are shown in Fig. 3 together with the calculated pattern for the 0-tt JJs. 

The magnetic field dependencies I C ,\(H) and I Cj2 (H) of the reference junctions are 
nearly ideal Fraunhofer patterns with oscillation period HqAH 93 //T, see Fig. 3. The 
small shifts along the if -axis are attributed to a net magnetization of the F-layer. For 
the 0-7T JJs the I C (H) oscillation period was ^184 /iT, nearly twice as large as for the 
reference JJs, as expected from theory. The plateau with a weakly developed dip in 
Fig. 3(a) and partially developed dip in Fig. 3(b) at H « are caused by the critical 
current asymmetry, i.e., l c \ 7^ —I c ,2- This can be seen from the calculated I C (H) pattern, 
too. For both samples the asymmetry of the main maxima and the shift of the I C (H) 
pattern along the if -axis indicate a difference in the net magnetizations $m,«- The best 
fit of the experimentally measured curves using Eq. (1) with $m,; as fitting parameters 
yields the fluxes $m,« < 0.2$ 0; i- e - ; dividing them by an area 6nm x 50 /mi of the F- 
layer in one half of the junction we obtain the net magnetization < 1.4 mT, whereas 
for a fully polarized NiCu alloys the saturation magnetization of 100 mT was reported 
[19]. The fitting parameter gives the mean magnetization in each half, indicating that 
the zero field cooled F-layer is in a multi-domain state. We attribute the remaining 
discrepancy with data to local nonuniformity of magnetization. Another feature - 
bumped minima in I C (H) for J C) i 7^ — i Ci2 and $m,i 7^ $m,2 [13] — appear very close to 
the measurement resolution limit and have not been further investigated. The I C (H) 
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dependence of the sample with c? 2 = 8.13nm shown in Fig. 3(c) is very symmetric, 
indicating $m,i = $m,2 = 0. The maximum critical currents at the left and right 
maxima of the I C {H) pattern for the 0-tt JJ are 99.2 and 98.9 fiA, respectively. They 
differ by less than 1%, and were ~ 0.72I ci , as expected from the theory [14]. The 
central feature of 0-7T JJs — the dip or plateau at zero field with critical current given 
by \I c (0)/I Ct i\ = (1 — \I C ,2/Ic,i\) /2 — is fairly affirmed for all three 0-7r JJs. 

In summary, we have fabricated SIFS JJs with and without a step in the F-layer 
thickness using state-of-the art SIFS JJ technology and produced 0-0, 0-tt, n-n JJ with 
different I c in the two parts. The experimentally measured I C {H) dependences can be 
well described by the model assuming different j c and different net magnetization in 
different parts of F-layer. The fits are much better than for the previous generation of 
samples [13]. The presented results demonstrate a good understanding of the physics of 
such 0-7r JJs. 

Thus the presented SIFS technology may provide the JJs with step- wise tailored j c {x) 

including the tt regions with j c (x) < 0. 
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Figure 1. (Color online) (a) The dependence \I c {<1f)\ ((<^f) as deposited initially) for 
JJs with non-etched (stars) and etched (circles) F-layer. The symmetric 0-ir junction 
(solid line) has I c (di) = ~I c {d2). (b) Sketch of a step-type 0-ir J J and its local phase 
<^i,2(a0- H an d ^M,i are orientated along the same axis. 




Figure 2. (Color online) Surface plot of I C {H) for I c 2/^,1 = — 1 • • • 1 and fluxes 
&M,i = calculated by Eq. 1. The measured I C (H) pattern (right scale, dotted line: 
baseline) were shifted by their I c .2/I c ,i ratio denoted in the black box (left scale). 
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Figure 3. (Color online) Measured I C (H) (voltage criteria of 0.5 /xV) of 0, ir and Q-w 
JJs for three J J sets. In (a) and (b) the 0-tt JJs have I Ct \ ^ —I c ,2 (dip at H « is not 
fully developed) . The asymmetry of the main maxima indicates a difference in local 
magnetizations 3>M,i- In (c) the 0-7r JJ is symmetric, i.e. I c _\ = —I c ,2- The insets 
depict the central dips of the asymmetric 0-7r JJs. Calculations (a)-(c) were done with 
&M,i = (dashed lines) and fitted $m,?: (solid lines). 



